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Abstract 

We present a preliminary measurement of the time-dependent CP- violating asymmetry in 
the neutral B decay B° — ► with — > iC 1- K~ and — ► 7r + 7r~. The measurement 
uses a data sample of about 87 million T(4S*) — > BB decays collected between 1999 and 
2002 with the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. In this 
sample we study events in which the CP final state is fully reconstructed and the flavor of 
the other neutral B meson is determined from its decay products. The amplitude of the 
CP-violating asymmetry sin2/3 is derived from the decay-time distributions. We measure 
sin2/3 = — 0.19±g;|o(stat) ± 0.09(syst). 
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1 Introduction 

Recent measurements of the CP-violating asymmetry parameter sin2/3 by the BABAR [jl|] and Belle [||] 
collaborators established CP violation in the B° system. These measurements, as well as the updated 
measurement of sin 2(3 = 0.741 ± 0.067(stat) ± 0.033(syst) by BABAR [^] reported at this conference, 
are consistent with the Standard Model expectation based on measurements and theoretical estimates of the 
Cabbibo-Kobayashi-Maskawa quark-mixing matrix [Q]. 

Charmless hadronic B meson decays provide important information for the study of CP violation ef- 
fects. The charmless B meson decays into final states with a meson are interesting because they are 
dominated by b — > sss gluonic penguins (Figure [I]), with a smaller contribution from electroweak penguins, 
while other Standard Model contributions are highly suppressed. These decays allow the extraction of the 
CP-violating parameter sin2/3. Comparison of the value of sin2/3 obtained from these modes with that from 
charmonium modes probe for new physics participating in penguin loops []|, ^]. The predicted deviation of 
the effective sin2/3 for the (f)Kg mode from sin2/? in the Standard Model is smaller than 4% [^,0]. In this 
analysis we probe for sizable deviations which are possible in many scenarios beyond the Standard Model. 




u, d u, d u, d u, d 



(a) (b) 



Figure 1: Quark-level diagrams describing the decays B — > cf)K: (a) internal penguin, (b) flavor-singlet 
penguin. 

The decay of neutral B mesons to the CP = - 1 final state (j)K® has been observed by BABAR in a sample 
of about 45 million B mesons with a branching fraction of BF(B° -> 4>K°) = (8.ll|g ± 0.8) x 10" 6 [§]. 
The channel B + — * 4>K + , which is used as control channel for the time-dependent analysis, was also 
observed with a branching fraction of {7-7t\l ± 0.8) x 10~ 6 [@]. 

The measurement of the time-dependent CP asymmetry in B° — ► (ftK® is similar to our approach in the 
charmonium channels We use an extended parametrization of the likelihood describing the event yield 
in signal and background which is combined with the likelihood for the decay-time distributions. 



2 The BABAR Detector and Data Set 



This measurement is based on data recorded with the BABAR detector [10] at the PEP-II energy-asymmetric 
e + e~ storage ring at SLAC. The data sample corresponds to an integrated luminosity of approximately 
80 fbr 1 that was collected at the T(4S) resonance. 

The detector consists of a five-layer silicon vertex tracker (SVT), a 40-layer drift chamber (DCH), a 
detector of internally reflected Cherenkov light (DIRC), an electromagnetic calorimeter (EMC), assembled 
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from 6580 CsI(Tl) crystals, all embedded in a solenoidal magnetic field of 1.5 T and surrounded by an 
instrumented flux return (IFR). The performance of the detector is discussed in [^]. 



3 Analysis Method 

3.1 Time dependent Analysis 

Each candidate event consists of a fully reconstructed neutral B meson, Bqp, decaying into (f>K® and a 
partially reconstructed recoil B, B tag , which we examine for evidence that it decayed as B° or B° (flavor 
tag). The decay-time distribution of B decays to a CP eigenstate with a B° or B° tag can be expressed in 
terms of a complex parameter A that depends on both the B°-B° oscillation amplitude and the amplitudes 
describing B° and B° decays to this final state [11]. The decay rate f+(f-) when the tagging meson is a 
B°(B°) is given by 



-\At\/r B o 

f±(Ai) = 



4t b o 



22mA . A . 1 - |A| 2 

1 ± - sin (Am d At) t : , ,,, 2 cos (Am d At) 

1 + | A| ^ l + |Ar 



(1) 



where At = tcp — ttag is the difference between the proper decay time of the reconstructed B meson (Bcp) 
and the proper decay time of the tagging B meson (B tag ), t b o is the B° lifetime, and Am^ is the B°-B° 
oscillation frequency. The sine term in Eq. ([j]) is due to the interference between direct decay and decay 
after flavor change, and the cosine term is due to the interference between two or more decay amplitudes 
with different weak phases. Evidence for CP violation can be observed as a difference between the At 
distributions of B°- and P°-tagged events or as an asymmetry with respect to At = for either flavor tag. 

In the Standard Model and for the case that the decay proceeds purely via b — > sss gluonic penguin 
transitions, A = r]fe~ 2t ^, and the angle (3 of the Unitarity Triangle of the three-generation CKM matrix [12] 
is given as (3 = arg [-V cd V* b /V td V t l\. 

Thus, the time-dependent CP-violating asymmetry is 

Acp{At) = f+(Atj +f"(At) = -7 ?/ sm2 /^ sm ( Ara rf At) , (2) 

with CP eigenvalue rjf = —1 for <j>K®, and K® —* 7r + 7r~. 
3.2 Reconstruction of the 4>K ° Final State 

We fully reconstruct B meson candidates (Bcp) in the decay mode 4>K® with Kg — > tt + tt~ and <fi — > 
K + K~ . For the charged tracks belonging to the recoil B and to the we require at least 12 measured 
drift chamber hits and a minimum transverse momentum of 0. 1 GeV/c. The kaon tracks of the <p in addition 
have to originate within 1.5 cm in the transverse plane and 10 cm in beam direction from the interaction 
point. 

A track is identified as a kaon based on a likelihood ratio combining the dE/dx information from 
the SVT and DCH below 700 MeV/c, and from DCH dE/dx and DIRC Cherenkov angle and measured 
Cherenkov photon number above this momentum. We define <fi candidates as pairs of tracks with opposite 
charge which can be combined and fit to a common vertex and whose invariant K + K~ mass lies within a 
20 MeV/c 2 mass interval centered at the <ft mass. In case there is no particle identification (PID) information, 
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the kaon hypothesis is assumed for one of the two tracks from the candidate (f> decay. Using a relativistic 



Breit-Wigner function of fixed mass and width [ |13| ] convoluted with a Gaussian we obtain an invariant mass 
resolution of 1.1 MeV/c 2 in the (f> signal. 

Analogously to the <p, we construct the from two oppositely charged tracks which are assumed to 
be pions. The selection of Kg — > tt + it~ candidates is based on the angle a between the line connecting 
4> vertex and K^-decay vertex and the momentum direction reconstructed from the pions (cos a > 0.999). 
Furthermore, we use the decay-time significance t/a t (t/<?t > 3). 

3.3 Event Yield Variables 

The measurement of B decays at the T(45) resonance provides kinematic constraints for the initial state. 
Substitution of the measured energy by the beam energy reduces the resolution of kinematic variables sub- 
stantially. 

Energy resolution can be expressed as: 

AE = E B - E bc , (3) 

with E\j C the beam constrained energy, which for the candidate B meson is derived as follows: 

_ S + 2pj • p B 

Ebc - 2Ei ' (4) 

with yfs the total e + e~ center-of-mass energy. The four momentum of the initial state is represented by 
(Ei,pi), and (Eb,Pb) is the four momentum of the candidate B meson, both measured in the laboratory; 
Eb c results from the assumption that we have particle-antiparticle production. Notice that the l?-candidate 
momentum ps is independent of the mass values assigned to the tracks comprising the candidate B. Signal 
events distribute in AE according to a Gaussian with a mean consistent with zero (±2 MeV). The observed 
width is about 17 MeV. The background shape in AE is parametrized by a linear function. We require 
\AE\ < 200 MeV. 

The second kinematic quantity in our analysis is the beam-energy substituted mass ities, which is 
defined as: 

m ES = \JeIc~Pb- (5) 

Signal events are distributed Gaussian-like in ities with a mean at the B mass and a resolution of about 
2.6 MeV/c 2 , dominated by the beam energy spread. The background shape in ities is parametrized by a 



threshold (ARGUS) function [14] with a fixed endpoint given by the average beam energy. Our selection 
requires ities > 5.22 GeV/c 2 . 

The helicity angle Oh of the is defined as the angle between the direction of the decay K + and 
the parent B direction in the <p rest frame. For pseudoscalar-vector B decay modes, angular momentum 
conservation results in a cos 2 Or distribution. In this variable the background is uniformly distributed. 

Monte Carlo simulation demonstrates that contamination from other B decays is negligible. Possible 



(K + K~) S-wave contributions (/o(980)) are not expected to peak under the <p meson [15] and are sup- 
pressed by the helicity angle which distributes uniformly for this background. However, charmless hadronic 
modes suffer from backgrounds due to random combinations of tracks produced in the quark-antiquark (qq) 
continuum, where q is dominantly u, d, and s quarks. The distinguishing feature of such backgrounds is 
their characteristic event shape resulting from the two-jet production mechanism. 

We consider the angle 9t between the thrust axis of the B candidate and the thrust axis of the rest of 
the event, where the thrust axis is defined as the axis that maximizes the sum of the magnitudes of the 
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longitudinal momenta in the T(4S) center-of-mass system. This angle is small for continuum events where 
the i?-candidate daugthers come from back-to-back qq jets, and is uniformly distributed for true BB events. 
In the event preselection we require | cos 9t\ < 0.9. 

Additional shape information comes from the momentum flow around the B thrust axis described 
through momentum weighted Legendre polynomials, Lj. The best separation between the signal and contin- 
uum events is achieved with the zeroth order (Lq = P*) and the second order (L2 = J2Pi x ^( cos2 ®t~ 1)) 
polynomials, where p* and 9\ are the center-of-mass momentum and angle with respect to the Bcp thrust 
axis and the sum is over all charged tracks and neutrals in the event that are not associated with the B- 
candidate. 

The last event shape variable is the B production angle, 9p, with respect to the beam direction in the 
T(AS) center-of-mass frame. In decays of a real T(45) into two pseudoscalar B mesons, the production an- 
gle follows a sin 2 9b distribution, while it is approximately uniformly distributed for the continuum events. 

The shape variables are strongly correlated and cannot be used independently in the likelihood calcula- 
tion. A Fisher discriminant is formed as a linear combination of the shape variables x = | cos 9t\, cos 9b, 
Lq, and L 2 : 

F = ^2liXi, (6) 

where coefficients 7^ are chosen such to make the maximum separation between the signal and continuum 
event distributions. The coefficients are calculated using Monte Carlo signal events and background events 
from data sidebands (0.1 < AE < 0.3). For the resulting signal Fisher distribution we use a bifurcated 
Gaussian distribution. The background Fisher shape is described by a sum of two Gaussian distributions. 

4 Tagging and Vertexing 

We use a B-tagging algorithm based on multivariate techniques to determine the flavor of B tag The al- 
gorithm relies on the correlation between the flavor of the b quark and the charge of the remaining tracks in 
the event after removal of the Bqp candidate. Separate neural networks are trained to identify primary lep- 
tons from semileptonic B decay, kaons, soft pions from D* decay, and high-momentum charged particles. 
The outputs of each neural network are combined to produce five hierachical and mutually exclusive tag- 
ging categories. Events with an identified electron or muon, and a supporting kaon, if present are assigned 
to the Lepton category. Events with an identified kaon and a soft-pion candidate with opposite charges are 
assigned to the Kaon I category. Events with one or more kaon candidates, and no lepton or soft-pion can- 
didates, are assigned to the Kaon I or Kaon II categories depending on the estimated mistag probability. 
Events with only a soft-pion candidate are assigned to the Kaon 1 1 category as well. The remaining events 
are assigned to the Inclusive or Untagged category based on estimated mistag probability. 

The quality of tagging is expressed in terms of the effective efficiency Q = J2c e c(^ — 2w c ) 2 , where 
e c and w c are the efficiency and mistag probability, respectively, for events tagged in category c. Ta- 
ble [I] summarizes the tagging performance in a data sample of fully reconstructed neutral B decays into 
D ( -*>-h + (h + = 7T+, p+, af) and J/ipK*° (K*° -> A' + vr~) flavor eigenstates (B flav sample). The recoil 
B, B tag , is again partially reconstructed. We use the same tagging efficiencies and dilutions for the (ftK® 
channel extracted from the statistical dominant flavor sample. 

The time difference At is obtained from the measured distance between the z positions of the Bqp 
and B tag decay vertices and the known boost of the e + e~ system. For the Bqp we achieve a z- vertex 
position resolution of better than 60 \im which compares well to the resolution obtained in the final state 
J/^fKg. The z position of the B tag vertex is determined with an iterative procedure that removes tracks 
with a large contribution to the total % 2 . An additional constraint is constructed from the three-momentum 
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Table 1: Tagging efficiency e, average mistag fraction w, mistag difference Aw = w(B°) — 
w(B ' ), and effective tagging efficiency Q for signal events in each tagging category. The 
values are measured in the B fi av sample. 



Category 






w 


;%) 


Aw( 


%) 


Q( 


%) 


Lepton 


9.1 


±0.2 


3.3 


± 


0.6 


-1.5 ± 


1.1 


7.9 ± 


0.3 


Kaon I 


16.7 


±0.2 


10.0 


± 


0.7 


-1.3 ± 


1.1 


10.7 ± 


0.4 


Kaon II 


19.8 


±0.3 


20.9 


± 


0.8 


-4.4 ± 


1.2 


6.7 ± 


0.4 


Inclusive 


20.0 


±0.3 


31.5 


± 


0.9 


-2.4 ± 


1.3 


2.7 ± 


0.3 


Untagged 


34.4 


±0.5 

















Total Q 28.1 ± 0.7 



and vertex position of the Bqp candidate, and the average e + e~ interaction point and boost. For 98% of 
candidates with a reconstructed vertex the r.m.s. Az resolution is 180 fim (1.1 ps). We require \ At\ < 20 ps 
and <7 At < 3.5 ps, where a At is the event-by-event error on At. 

The At resolution is dominated by the tag-side (B tag ), which is well under control from our flavor 
sample. The empirical At resolution function for signal candidates is parametrized by a sum of three 
Gaussian distributions (see Ref. JI^]), with parameters determined from a fit to the flavor sample. A common 
parametrization is used for all tagging categories, and the parameters are determined simultaneously with 
the CP parameters in the maximum likelihood fit. The tagging parameters for the untagged events are fixed 
to w = 0.5 and Aw = 0. The At background in the (j)K® channel does not show a lifetime component and 
is parametrized as a sum of two Gaussians, core and tail, with the tail fraction of 2%. 

The parametrization of At is checked by measuring the lifetime in the channel <j)K + ( 1 80 signal events) 
which has a compatible At distribution. The obtained value agrees within la with the world average [|I3|]. 



4.1 Maximum Likelihood Fit 



We use an unbinned extended maximum likelihood fit to extract yields and CP parameters from the 4>K® 
(CP) sample, and simultaneously At resolution parameters and tagging quantities in the flavor sample. The 
likelihood for candidate j tagged in category c is obtained by summing the product of event yield rtj, tagging 
efficiency ej >c , and probability Vi )C over the two possible signal and background hypotheses i. 



AM eXP 



i 



n 



(7) 



The probabilities Vi yC are evaluated as the product of PDFs for each of the independent variables Xj = 
{mES; T, cos 6h, At} in the CP sample, and xj = {m^, AE, At} in the Bfi av sample. The c?j are 
fixed parameters that describe the expected distributions and are derived from fits to signal Monte Carlo, 
the B + — > 4>K + control channel, on-resonance sidebands, and off-resonance data. The At distributions of 
the Bflav sample evolve according to flavor oscillation in B° mesons. The observed amplitudes for the CP 
asymmetry in the CP sample and for flavor oscillation in the flavor sample are reduced by the same factor 
1 — 2w due to flavor mistags. The total likelihood C is the product of likelihoods for each tagging category 
and the free parameters are determined by minimizing the quantity — In C [17]. 
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The total number of (f>Kg candidate events in the fit region is 1352. In order to extract the event yields 
we perform an initial fit without tagging or At information. Our final CP sample is composed of 5 1 signal 
and 1301 background events distributed over the fit range. For our CP fit we fix these yields. The Bfi av 
sample consists of about 26000 signal events with a purity of better than 80%. 

There are 34 free variables in the fit, in agreement with the charmonium sin 2(3 fit [§], where only the 
parameter sin 2(3 is solely fit from the signal in B° — > 4>K® (|A| = 1 fixed). The other parameters are the 
average mistag fraction w and the difference Aw between B° and B° mistags for each tagging category (8), 
parameters for the signal At resolution (8), parameters for background time dependence (6), background 
At resolution (3), and mistag fractions (8). The determination of the mistag fractions and At resolution 
function parameters for the signal is dominated by the high-statistics flavor sample. We fix t b o = 1.542 ps 
and Amd = 0.489 ps -1 [|l3|]. The largest correlation between sin2/3 and any linear combination of the other 
free parameters is 2%. 

The result for the effective sin 2/3 is: 

sin 2(3 = -0.19ig;§ (stat) ± 0.09(syst) (8) 
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Figure 2: The distribution of the (j)Kg events in four fit variables with tighter selection criteria corresponding 
to the ranges shown (see text). The tues distribution is presented in a wider range. The solid line refers to 
the fit for all events, the dashed line corresponds to the expected background distribution. 

Figure ^| shows the event yield variables for all events in the limited ranges 5.27 < rriEs < 5.3 GeV/c 2 , 
\AE\ < 0.1 GeV, — 2 < T < 3, and | cos##| > 0.2, focussing into the signal region. Figure || shows the 
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At distributions for the B° and the B° tagged subsets of these events with the fit superimposed. 

As a consistency check we measure the CP asymmetry in the final state B + — ► <fiK + . From a sample 
of 180 signal events we obtain a CP asymmetry of sin 2fi{4>K + ) = 0.26 db 0.27, which is consistent with 
the expected value of zero. 

B° tags: 



c 

> 




B° tags: 



4 6 
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d 



a 
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pq 



BaBar 

preliminary 




Figure 3: Time distributions for B and B° tags. The solid line refers to the fit for all events, the dashed line 
corresponds to background. 

Repeating the fit with all parameters except sin2/3 fixed to their values at the maximum likelihood, we 
attribute a total contribution in quadrature of 0.01 to the error on sin2/3 due to the combined statistical 
uncertainties in mistag rates, At resolution, and background parameters. 

The fit is repeated on generated datasets based on the probability density functions for signal and back- 
ground shapes. We do not observe a bias in the refit value of sin 2/3. 

We consider systematic uncertainties due to the event yield determination in 4>K® (0.02), limited Monte 
Carlo statistics (0.02), composition and CP asymmetry in the background in the CP events (0.03), the 
assumed parametrization of the At resolution function (0.02), due in part to residual uncertainties in the 
Silicon Vertex Tracker alignment, and the fixed values for Am^ and tb (0.006). 

Furthermore, we explore the sensitivity to the parameter |A| in our limited sample. It turns out that the 
fit is not simultaneously sensitive to both sin 2/3 and |A|. Therefore, we scan for the value of |A| over a 
wide range (| A| = 0...3) but do not observe a strong variation of the central value of sin 2/3. We attribute an 
additional conservative error of 0.08 to our value of sin 2/3, which is the maximum variation observed in the 
scan. 
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5 Conclusions 

We measure the preliminary effective value of the time-dependent CP asymmetry sin 2/3 = —0 .19~^Q^(stat) 
±0.09(syst) in the decay of neutral mesons into the final state 4>K®, K® — ► tt + it~. The deviation of 
this value from the updated BABAR value presented at this conference, sin 2/3 = 0.741 ± 0.067(stat) ± 
0.033(syst), is about two standard deviations. The measurement will, for some time, be dominated by the 
statistical uncertainty. 
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